Flow equation approach to heavy fermion systems 
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We use Wegner's flow equation method to investigate the inflnite-C/ periodic Anderson model. 
We show that this method poses a new approach to the description of heavy fermion behaviour. 
Within this scheme we derive an effective Hamiltonian in which the c and / degrees of freedom are 
decoupled. By evaluating one-particle energies as well as correlation functions we find an electronic 
structure which comprises two gapped quasiparticle bands. We also address the lattice Kondo 
temperature, which shows a typical exponential dependence on the hybridisation energy. This 
energy scale exhibits a significant decrease compared to that of the single impurity Anderson model. 

PACS numbers: 71.10.Fd,71.27.+a,75.30.Mb,71.20.Eh 



I. INTRODUCTION 

Heavy fermion systems (HFS) have triggered a grow- 
ing interest in both experimental and theoretical physics 
over the last decades. They are mainly based upon lan- 
thanide or actinide compounds, thus their physics is dom- 
inated by the interplay of itinerant (3d) electrons and 
rather localised electrons from partly filled 4/ or 5/ elec- 
tron shells. These systems generally posses rich phase 
diagrams, and they show intriguing physical phenom- 
ena, such as superconductivity, long-range magnetic or- 
der, Fermi-liquid and non-Fermi-liquid behaviour, and 
in many cases magnetism and superconductivity coex- 
ist. Despite much theoretical effort, many open ques- 
tions remain even for the paramagnetic metallic phase. 
At low temperatures this phase is generally characterised 
by the existence of heavy quasiparticles, which lead to a 
strongly enhanced specific heat coefficient as well as an 
enhanced susceptibility. Experiments suggest that the 
effective mass of the quasiparticles is in the region of 
several hundreds up to thousand bare electron masses. 
On a microscopic basis, HFS can be described by the 
periodic Anderson model (PAM) and the Kondo lattice 
model (KLM). Whereas the first describes a system of 
localised / electrons and conduction electrons, which in- 
teract via a hybridisation, the latter yields the micro- 
scopic picture of a periodic range of magnetic moments 
in a metallic host. In the so-called Kondo limit the KLM 
can be regarded as an effective model for the PAM. The 
physics of both models is associated with the Kondo ef- 
fect [l| , which was initially examined in connection with 
a single magnetic impurity interacting with a sea of free 
conduction electrons. The typical energy scale for this 
interaction is provided by the Kondo temperature Tk, 
which shows an exponential dependence on the hybridi- 
sation strength. This relation cannot be understood in 
terms of conventional perturbation theory. Beside sys- 
tems that involve / electrons, heavy fermion behaviour 
has also been observed in LiV204 [2|- It is believed that 
the occurrence of a heavy Fermi liquid state is driven by 
the frustration of the intersite spin couplings 3] . 

About 10 years ago Wegner |j| and independently 
Glazek and Wilson '91 proposed a scheme with the objec- 



tive to diagonalise or block-diagonalise Hamiltonians on 
the basis of a continuous unitary transformation. This 
method is known as the flow equation method and has 
been used for the theoretical description of a great va- 
riety of physical systems. In a recent application this 
method was used to investigate the Hubbard model, a 
paradigm for a strongly correlated lattice system Q . Al- 
though the authors treated the system on the basis of 
perturbation theory, they obtained excellent results in 
the regions of moderate and even strong coupling. In this 
paper we apply the flow equation method to the infinite- 
U PAM in order to describe the physics of the paramag- 
netic metallic phase of HFS. Thereby we do not consider 
perturbative arguments with respect to the hybridisation 
strength. We focus on the electronic structure, which 
suggests the existence of heavy quasiparticles, and show 
that the lattice Kondo temperature features an exponen- 
tial dependence on the hybridisation strength. This pa- 
per starts off with some general remarks regarding the 
method and the microscopic model. Afterwards we dis- 
cuss the application of the method as an approach to the 
PAM. The last part is devoted to a detailed discussion of 
our findings. 



II. MODEL AND METHOD 

A. Periodic Anderson Model 

The PAM is the standard model for the description of 
heavy fermion systems. Supposing the Coulomb repul- 
sion of the local / electrons being the dominant energy 
scale in the system, this model can be written in the form 



iifpAM = 51 ^k cl^c^^ + ^£f flu 
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In this representation fj^ creates an / electron on site i 
with spin a. Due to the infinitely strong Coulomb repul- 
sion of two / electrons doubly occupied sites are explicitly 
excluded. This is done by the use of Hubbard operators 



rather than usual fermionic creation and destruction op- 
erators: 



fi = /I n(i 
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(2) 
(3) 



Here n{^ denotes the number operator for the / elec- 
trons. The essential property of the Hubbard opera- 
tors is that the creation operator fj^ only acts on an 
empty site i, and the operator /^^ only annihilates a 
singly occupied state. The operator Pf stands for the 
projector on the empty / site i. In this context the 
physical subspace merely comprises singly occupied and 
empty /-sites. As a consequence the completeness rela- 
tion P? -I- Ylicr fiaha = -*- holds, and the operators obey 
the following anti-commutation relations 



\/i(T' h'a' /+ ^ "m' y^Ta' Jia-Jia + JiaJirr') 



(4) 



The Fourier transform of the Hubbard operators 
can be defined in the usual way as f^^ = 
1/viV^^ exp(ikRi)/2^, where N is the number of lat- 
tice sites. However, the above introduced commutation 
relations for the Hubbard operators in the local picture 
result in rather complex commutation relations for the 
operators f^^ and f^^. In Eq. Q cj^.^ creates a conduc- 
tion electron with momentum k and spin a. Whereas the 
/ electrons are localised, the conduction electrons possess 
a dispersion e^. The two subsystems are connected via 
a hybridisation of the strength V^. Although, in general, 
a dispersive hybridisation may be examined, we here fo- 
cus on the case Vk = const, i. e. an interaction which is 
restricted to electrons on the same lattice site only. As 
a further simplification we only consider a hybridisation 
between c and / electrons within the same spin channel. 
The model (^ has been in the focus of a great variety 
of theoretical investigations. One of the most descriptive 
ones is the slave-boson mean-field (SBMF) theory Q. 
Here the Hubbard operators are treated by the introduc- 
tion of auxiliary boson fields, which ensure the conserva- 
tion of the commutation relations. On the other hand, 
the configuration space of the particles must be restricted 
due to their bosonic character. This can be accomplished 
by considering the constraint with the help of additional 
couplings in the Hamiltonian. The used Lagrange mul- 
tipliers have to be determined self-consistently. For this 
model the problem can be solved in the limit i/y — > oo, 
where the latter quantity denotes the spin degeneracy of 
the / orbitals. A quite similar treatment of the model is 
provided by the Gutzwiller variational method ^ , which 
can also be formulated in terms of auxiliary bosons Q. 
The intriguing result of this approach is an instability 
of the ground state with respect to ferromagnetic order- 
ing for the case Vf = 2, and the Kondo temperature is 
significantly enhanced compared to the single impurity 



Though the analytical approaches to the PAM have 
given insight to various aspects of heavy fermion sys- 
tems, they generally have to deal with more or less rude 
simplifications. In the last decade a range of numerical 
solutions emerged. A great number of recent works are 
based on the Dynamical mean-field theory (DMFT) [lO| . 
Here the PAM can be mapped onto an effective single 
impurity Anderson model (SIAM) . This approach proves 
to be exact in infinite dimensions, and it is considered to 
deliver good results even for three-dimensional systems. 
The major task of this scheme is a self-consistent numer- 
ical solution of the impurity problem, whereby a number 
of methods, such as Wilson's numerical renormalisation 
group (NRG) [ij. Quantum Monte Carlo (QMC) simu- 
lations [l^ and iterated perturbations theory (ITP) i;||, 
have been used. Other numerical approaches comprise 
the diagonalisation of small clusters 14] and QMC [13 
for the full system. As the numerical results of Sec.lVlare 
obtained for a one-dimensional system, we should men- 
tion that this system and the corresponding KLM have 
been extensively studied in the framework of the Density 
Matrix Renormalisation Group (DMRG) 16]. Various 
studies of the KLM suggest this system being in the uni- 
versality class of Luttinger liquids. The phase diagram 
of the one-dimensional PAM was subject of a series of 
recent works 17]. 

Some basic properties of the model Eq. ^ have also 
been discussed in the framework of a renormalisation ap- 
proach by Becker et al. |18| , which incorporates the use of 
hard cut-off functions. However, the authors use pertur- 
bation theory with respect to the hybridisation strength 
and neglect processes that connect different / sites, which 
is not sufficient to obtain the proper heavy fermion be- 
haviour. 



B. Flow equation method 

The flow equation method is based on the formulation 
of a continuous unitary transformation. For this rea- 
son in various works this approach is refered to as CUT. 
The transformation is represented by a unitary Opera- 
tor U{£), for a given Hamiltonian H. Here £ denotes 
the flow parameter 19] and the according transformed 
Hamiltonian is simply given by H{i) — U{t)HU'^{tj. By 
definition the operator C/(l) stands for the identity oper- 
ator and the fully transformed Hamiltonian H{(. = 00) is 
written as H* . The latter may be regarded as an effective 
Hamiltonian, and due to the unitarity of the mapping it 
possesses the same system of eigenenergies as H . The 
flow of H{t) can be conveniently formulated in its differ- 
ential representation 



dffffl 



[n{£),H{£)\ 



and it is governed by the generator 

dU{£) 



m 



A£ 



U\£) = -f]\l) 



(5) 



(6) 



This operator can be chosen to map H onto a much sim- 
pler Hamiltonian. If one starts out, for instance, with a 
Hamiltonian of the structure H = Hq + Hi, where Hq 
represents a system of free particles and Hi an interac- 
tion, the choice 



ri=[Hf,,Hi] 



(7) 



indeed triggers the interaction part to vanish in the limit 
i' -^ oo. This can be proved by some simple calculation 
and Eq. Q represents the original choice of the gener- 
ator by Wegner [J|. In this context the resulting effec- 
tive Hamiltonian is diagonal or block-diagonal. Beside 
Wegner's generator, in the recent years a variety of di- 
verse generators has been used. As we have mentioned 
in the introductory section, an independent approach for- 
mulated by Glazek and Wilson was made, in which the 
Hamiltonian is diagonalised by integrating out higher en- 
ergy contributions. Another choice for the generator was 
proposed by Knetter and Uhrig _2Qj. It has the intrigu- 
ing feature that it allows the construction of a particle- 
number conserving effective Hamiltonian. 

Independent of the choice of any particular genera- 
tor, the transformation can as well be applied to formu- 
late the flow of observables. Similar to the structure we 
have derived for the Hamiltonian, an arbitrary operator 
is transformed as 



dA{£) 

de 



HeiAie)] 



(8) 



As a result of this we can take advantage of the invari- 
ance of expectation values and correlation functions un- 
der unitary transformations. In this sense these quanti- 
ties can be evaluated with respect to the much simpler 
Hamiltonian H* by the simultaneous transformation of 
the operators. If we regard for example the expectation 
value of an observable A, the relation 



(^> 



TiAe-l^" T^^*g-/3/f* 



(9) 



Z Z 

holds, where Z is the partition function. 

The above considerations are instructive for the eval- 
uation of Green's functions. If we address, for instance, 
the retarded Green's function {{A{t); B{t'))) — —i9{t — 
t'){{A{t),B{t')}+), we can use the fact that these quan- 
tities are invariant with regard to a unitary transforma- 
tion. This allows us to evaluate the expectation value 
and the Heisenberg time dependence of the operators 
with respect to the transformed Hamiltonian H* , if we 
use the transformed operators A* and B* at the same 
time. In terms of the Fourier transformed Green's func- 
tion this can be simply denoted as {{A;B)){uj + iS) — 
*{{A*;B*))*{uj + iS). 



III. FLOW EQUATIONS FOR THE PERIODIC 
ANDERSON MODEL 



After having discussed the general aspects of the model 
and the method we now focus on the application of the 



flow equation approach to the infinite-/7 PAM. This sec- 
tion is organised as follows. First we introduce an ap- 
propriate generator with the objective to map the model 
onto a much simpler one, and subsequently we discuss 
the flow of the Hamiltonian. In a second step the trans- 
formation of observables is examined and the evaluation 
of expectation values is discussed. The last part of this 
section then is devoted to the solution of the resulting 
flow equations. 



A. Flow of the Hamiltonian 

As we have mentioned in the introductory remarks re- 
garding the flow equation method an appropriate choice 
of the generator is provided by iQ. The aim of our treat- 
ment of the PAM is the decoupling of the distinct types 
of electrons, namely the c electrons and the / electrons. 
They are connected by a hybridisation which is charac- 
terised by the energy Vk. If we identify this contribution 
with the interacting part of the Hamiltonian, we recog- 
nise that the remaining contribution is not diagonal in 
the usual sense. This is due to the representation of the 
/ electrons by Hubbard operators and results from the 
incorporation of their strong Coulomb interaction. Con- 
sequently, we use for the generator the more general form 

^w = E^kW(cL/u.-/Lo (10) 

kCT 

The coefficient rj^ has to be determined later. A similar 
definition has already been used for the flow equations 
analysis of related models [2^ |22| • 

A first remarkable result of the above choice for the 
generator is the generation an / electron hopping of the 
form 






(11) 



as a consequence of the evaluation of the commutator in 
Eq. (O . A detailed discussion of all contributions is given 
below. The dispersion of the / electrons is obtained by 
the Fourier transform of the hopping matrix element Uj 
with the exclusion of local processes 



^k(^) = ^E*»^-(^)^"^'^^'''"''^'^ 



(12) 



''#j 



With this additional interaction term the Hamiltonian of 
the flow results in the following form: 

HpAuii) = E £k(^) cLck. + E(^/(^) + ^k(^)) flAa 



k(T 



E^k(^)(cL/k. + /LO + ^(^) (13) 



ktr 



Besides the hopping of the / electrons an additional en- 
ergy shift E{£) was introduced. 



As a general convention for the upcoming consider- 
ations we suppress the exphcit declaration of the (.- 
dependence of all energies. We only use it in some cases 
in order to avoid confusions. Furthermore, we use for 
a certain quantity Q the representation Q' for its initial 
value, i. e. Q{i. — 0), and Q* for the limit i — > oo. 

Eq. O poses the central relation of the flow equation 
method, thus our task is the evaluation of the commu- 
tator of the generator H1U|) and the Hamiltonian (|13|l . 
As the scheme suggests the resulting contributions de- 
termine the flow of the Hamiltonian. The commutator 
can be easily calculated by making use of the commuta- 
tor relation Q) for the Hubbard operators. By doing so 
we arrive at the following result: 



context the fourth line of H14|l is evaluated as 



l^^e^(k'-k)R.^^^^,^t^^^^^^o + l^. 



N 



kk' i(7 



2 E 77k^k(P^) cLck. 'nT.T. '^k^k^J Ilh 



k(T 



ik crcr' 



^E''k^k(nL)(i-(^"» 



N 



(15) 



zkfT 



h ^] = - E '^k (^k - e/) (cL/k. + /Lck.) 



k(T 



4^EE'?ki^,e-'^^*P°cU. + h-C. 



ko- i^j 



+ -7W E E'^kt., e-'''^-/i,4cL4, + h. c. 



kcrcr' i^j 



^EE^^^'''"''"''^>^^''-Uv^°+h. 



kk' ia 



^ E E ^^^'"'^''' ^^v^' 'L-^.JIL + h. c. 



kk' icra-' 



-^EE2e-*^"'-"^^'?'^^«^/l4 



(14) 



kfj- ij 



In the first term we recognise a contribution, which pos- 
sesses the form of the hybridisation interaction. The last 
term represents the generated hopping of the correlated 
/ electrons, whose generation we have already discussed 
in the last paragraph. In terms of Hubbard operators 
both contributions can be characterised as one-particle 
contributions, whereas the remaining part is comprised 
of two-particle interactions. The second and the third 
term both correspond to a hybridisation-like interaction, 
which couples to an empty and a singly occupied / site, 
respectively. Alike the usual hybridisation term they 
alter the number of particles in the respective electron 
subsystem by one. In contrast the interactions in the 
fourth and fifth line do not change the number of par- 
ticles within the conduction band or within the / elec- 
tron subsystem. The penultimate term can be regarded 
a Coqblin-Schrieffer interaction. 

To cope with all contributions of Eq. (|14|l we are 
obliged to introduce a truncation scheme. The objec- 
tive of this is to keep the Hamiltonian in the form p3|l . 
This task may be accomplished by various techniques, 
one of them is the normal ordering of the operators [J . 
Because normal ordering must be defined with respect to 
a bilinear Hamiltonian, it is more convenient to use an 
alternative scheme here. We simply decouple the higher 
interactions by use of a Hartree-Fock decoupling. In this 



where the representation (n^^) = {cj^^c^^^) and the rela- 
tion P^ + J2a ficrfia = -*- were used. Here the expectation 
values are to be evaluated with regard to the Hamilto- 
nian (Q, thus they are constant during the entire flow. 
Obviously, the first two terms on the right hand side 
contribute to the one-particle terms of the Hamiltonian, 
whereas the last one results in an energy shift. It should 
be emphasised that this factorisation scheme preserves 
the character of the initial interaction with regard to 
the number of particles in each electron subsystem. In 
Eq. p5|l . for example, both the interaction on the left 
hand side and the resulting terms on the right hand side 
preserve the number of c and / electrons. Furthermore, 
the local projector P^ remains untouched by such a kind 
of approximation. A detailed derivation of the decou- 
pling of the remaining interaction is given in Appendix 

HI 



If we collect all the terms that are provided by the 
factorisation scheme, we recognise that the occurrence of 
the expectation value (P°) is always accompanied by the 
expectation value (ni^). As we have mentioned above the 
first denotes the average number of empty / sites, and the 
latter represents the number of singly occupied / sites. 
Since we consider the non-magnetic phase in the present 
paper, we can deduce the relation (nj^) = (1 — {Pf))/i'f 
from the completeness relation for the local / configura- 
tion space. As a consequence of this we conclude that 
{Pf) is of order i/9 whereas (f^fg.) is of order z/7^. The 
same becomes obvious for the contributions for the en- 
ergy Ey , where the terms of order v'i come from the parts 

ex Pf , and the terms of order vl are given by the contri- 
butions oc fl^f^^. As a simplification we consider in the 
following the case of great degeneracy of the electrons, i.e. 
I'/ — > cxD, and neglect all terms of order vj . The spirit of 
this approximation is rather similar to the SBMF theory, 
though both approaches cannot be simply compared on 
this level. Eventually, we arrive at the following system 
of differential equations for the one-particle energies of 



the Hamiltonian (|13() 
deic 



dAk 

dg/ 
d^ 



-- - 277k ^k 



(16) 
(17) 



2 



TV- 



E«.)'?kl4 - ^ E(^k.)r/kAk (18) 

(19) 



k(T 



^=-(-k-g/-(i^")Ak),k 



l^(Ak.)(l-(/^'^))77kAk 



(20) 



iWa 



The system possesses the initial conditions ek(^ = 0) = 
4, e/(^ = 0) = 4, Vk(^ = 0) = T/i^ and Ak(^ = 
0) = -E(^ = 0) = 0. From the above equations we 
can immediately conclude (P°)Ak(^) = Ek — £k(^) and 
^(^) = N{e'j - e/(^))(l - (if)) so that Ak can be elim- 
inated completely, by making use of the first relation. 

At this point we emphasise that we have not used per- 
turbational arguments with respect to the hybridisation 
strength nor have we used a mean-field decoupling of 
the Coulomb interaction within the above approximation 
scheme. This fact is very important with regard to the 
description of heavy fermion behaviour. 

With the form-invariant Hamiltonian H13(l we know 
that in the limit i* — j oo the resulting effective Hamil- 
tonian is of the form 



H, 



PAM 



E* t 
^k^kcrCko 

k(T 



ko- 



K)fLf^. + E* 



(21) 



It describes two decoupled subsystems of uncorrelated 
and correlated electrons and can therefore not be re- 
garded as diagonal. This fact is the main difference to 
other approaches, such as the SBMF theory, which leave 
a noninteracting Hamiltonian. Another important differ- 
ence is that we do not need to incorporate a side condi- 
tion, which restricts the number of particles on each / 
sites, for we did not made use of auxiliary boson fields. 



B. Flow of Observables 

In the introduction we have already discussed the rep- 
resentation of expectation values within the framework 
of flow equations. There we concluded that expectation 
values as well as correlation functions can be easily eval- 
uated on the basis of the derived effective Hamiltonian, 
if we regard the transformed operators at the same time. 
The transformation of an arbitrary operator works quite 
similar to the transformation of the Hamiltonian. If we 
consider, in the first instance, the flow of the operator 



Ckcr, we can write the mapping in the form 



d£ 



[11, Cka\ 



(22) 



Similar to the transformation of the Hamiltonian we have 
to regard all the contributions that are generated by the 
commutator on the right hand side. The operator, of 
course, has to fulfil the boundary condition Ckcr(O) = Cko-- 
Hence we obtain the lowest order contributions simply by 
evaluating the commutator with Cko-- With this results 
we can infer the following form for the ^-dependent op- 
erator 



Cko-(^) = Ok Ckcr + /3k fka 



(23) 



This is a linear combination of the ordinary operators Cko- 
and /ko- with the weighting factors ak and /3k, respec- 
tively. It is the simplest possible representation for the 
flow of the operator Cko-, and higher interactions that are 
generated by the commutator must be truncated in the 
same manner as we have treated the flow of the Hamilto- 
nian I23 . The complete evaluation of the flow equation 
yields differential equations for the weighting factors 



dttk r. /r,o\ d/?k ,„.^ 

-^ = ~l3kVk{Pi) , — = ak?7k (24) 

Here again the ^-dependence has been omitted for the 
reason of readability. Obviously, the initial values of the 
weight factors are given by aj^ = 1 and P^. = 0, and as a 
consequence of unitarity they obey the relation 



al 



(^°>/3, 



1 



(25) 



In a further step the coefhcient /?k can be eliminated, and 
the flow of ttk is governed by the equation 



^arcsin[l-2a2]=2y^77k 



(26) 



The transformation of the operator /ko- can be treated 
likewise, and we obtain 



kAi)^-{P^)/3kCka+ak^ 



ko- 



(27) 



With the results of the previous paragraph we are now 
in the position to evaluate all expectation values that oc- 
cur in the differential equations by using the unitarity 
of the transformation, i. e. (rikcr) — *(cJ^|jCJct)*. Here 
*(•••)* refers to the expectation value with respect to 
the fully transformed Hamiltonian pif) . By inserting the 
transformed operators, we obtain for the non-local ex- 
pectation values 



("-kcr) 



("k)^*(cLckJ 



m'^ifLLr (28) 



and 



(Ak.) = 2alPl (-(/^°)*(cLck.)* + *(/L/. 



(29) 



In view of the Hamiltonian (|21|l it is clear that there 
are no contributions from expectation values of the form 

The derivation of the non-local expectation values, 
which contain contributions from different lattice sites, 
could be carried out rather intuitively, since we made 
use of a descriptive transformation scheme for the single 
operators. In contrast, the treatment of the local expec- 
tation value (P°) turns out to be more difficult. Though 
the projector can be represented in terms of Hubbard 
operators as P° = fJa-fi^, its transformation cannot sim- 
ply be carried out by the transformation of the single 
operators, for this procedure would break the restriction 
for the / electrons to empty and singly occupied sites. 
Instead of that we must transform the entire operator 
P°. The detailed derivation of this expectation value is 
presented in Appendix^ and as result we obtain 



kCT L 



^kcr^kcr/ 



+ (1 -(/?:)') *(/L/k.> 



(30) 



/ Green's function is simply given by 

^((/ic.;/L»(^ + *'5)=({/k.,/L}+> 

(32) 

where we set fi, = 1. This representation suggests a cou- 
pling to higher order Green's functions, which can be 
seen by the second term on the right hand side. In order 
to obtain a closed equation of motion we again apply the 
decoupling scheme of the last section and neglect contri- 
butions of order i^J^ ■ As result we obtain 



{{L-JL))i^ + ^s) = 



(pp 



u; + ts-{e} + {P^)Al) 



(33) 



and, by the application of the spectral theorem, the de- 
sired expectation value reads 

* {fLLr = {P") ^F {e} + (P°) A* ) (34) 

On the right hand side the Fermi distribution np — (1 + 
exp{(3E))~^ has been introduced. 



Altogether, this more sophisticated treatment can also 
be applied in order to derive the non-local expectation 
values, if we transform them as pairs rather than as sin- 
gle operators. The result of this evaluation, though, co- 
incides with the above derivation, and we can use the 
much easier representation from the last paragraph first 
for reasons of clearness and second as a starting point for 
the evaluation of the Green's functions. 

Further on, we wish to discuss the form of the expecta- 
tion value which is provided by (|3()|l . If we alternatively 
start out from the fact that the total number of particles 
must be conserved under the unitary transformation, we 
can derive this result in another way. The relation 



E< 



^L)+"E(/^4 



ko- 



(cLcj^^) 



Z^ Uko-Jkcr/ 
ko- 



(31) 
denotes the particle number conservation. If we use on 
the left hand side the expression H28|) for the expectation 
value (n^^) and divide by the number of sites, we eventu- 
ally arrive at Eq: (|30|) . Consequently, we can deduce that 
the above derivation is consistent with the conservation 
of the total number of particles. 

As a remaining task we have to evaluate the expec- 



ts ^ktr/ 



and *(/L/ko)* ■^it'^ respect to 



tation values *(c^ 

the effective Hamiltonian (I21|l . As the subsystem of c 
electrons simply consists of non-interacting fermions the 
respective expression is determined by a Fermi distri- 
bution. The calculation of the remaining /-subsystem 
is more complex as these particles do not obey a sim- 
ple statistics. It is expedient to calculate the demanded 
expectation value by means of ordinary Green's func- 
tion calculations |24l|. The equation of motion for the 



IV. INTEGRATION OF THE FLOW 
EQUATIONS 

A. Analytical solution 

The PAM during the flow is given by Eq. H13|l and leads 
to an effective Hamiltonian of the form (|21|l . provided 
we have chosen a generator which ensures the vanishing 
of the hybridisation Vu- The evolution of the matrix 
elements is given by a closed set of coupled differential 
equations, which have to be evaluated for the derivation 
of physical quantities. In a first step we try to find an 
analytical solution. For this purpose we bring Eq. 119|) 
in the following form 
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Here we have substituted A^ following the above dis- 
cussion. This equation provides a representation of 77k 
as function of the one-particle energies and allows us to 
eliminate this quantity in all other differential equations. 
With Eq. H35|l we are now in the position to rewrite the 
differential equation for the single-particle energy et as 
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As a simplification for the next discussions we assume 
that the energy ej converges rather fast to its limit value, 
while the main contribution to the evolution of the other 
energies occur at a larger £-scale. This assumption al- 
lows us to replace e/(^) by the constant energy el. Such 
kind of simplification is quite common in terms of flow 



equations |21ll22| . With this approximation the equation 
(|36|l can simply be integrated and we arrive at 

(elf - el{e} + 4) + 4^} - (^- >(K)' = (37) 

using the fact that Vk vanishes in the hmit ^ — > cx). This 
quadratic equation has the sohitions 



where we have used the representation 
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Equivalently, we obtain the foUowing solution for the / 
degrees of freedom 
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The latter describes the excitation energies of the / elec- 
trons as suggested by Eq. lIH^ . In order to derive a sim- 
ple equation for the energy £/ we neglect the term with 
(Ako-) on the right hand side of H18|l . and after having 
integrated the resulting equation, we obtain 
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(41) 
The simplification of a constant energy Sf can also be 
made for the derivation of the weighting factors and allow 
us the evaluation of expectation values. From (|26|l we 
arrive at 
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(42) 



Altogether the results coincide with the expressions from 
the SBMF theory. The excitation energies ||3HI) and (gOj) 
constitute the electronic structure of the system of two 
gapped bands. The overall renormalisation of the / level 
is provided by Eq. H41|l . These results are based on the 
assumption of a fast convergence of £/ compared to all 
other energies and give us a rough estimate of the elec- 
tronic structure of the full problem. However, a thorough 
examination beyond this simplification requires a numer- 
ical treatment of the differential equations. We would 
like to stress that we have not used a particular choice 
for the coefficient ?7k. Thus the derived results are valid 
for any ?7k that yields a vanishing hybridisation for the 
fully transformed Hamiltonian. 



B. Numerical solution 

The above derivation of an analytical solution for the 
differential equations was based on further simplifica- 
tions. A more general solution can be provided by a 
numerical integration. While the analytical solution was 



independent of a particular choice of the coefficient 77k, 
we now need to determine this quantity for the numerical 
evaluation. An expedient choice is given by 
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(43) 



since the energy Vk vanishes when t goes to infinity. This 
becomes obvious by inserting Eq. (|43|l into Eq. H19|) . for 
the differential of Vk is a negative quantity in the en- 
tire integration range, while Vk itself is always positive. 
The above choice is not dissimilar to Wegner's generator. 
However, it is not simply obtained by commuting the 
interacting and the non-interacting part of the Hamilto- 
nian, as suggested by Eq. (jT)) , in that the non- interacting 
part for the present calculations incorporates the hopping 
of strongly correlated objects. In fact, the derivation of 
the generator (|43|l uses a factorisation scheme, which has 
been presented in detail in Sec. IIII Al This particular 
form of the generator also triggers the Hamiltonian flow- 
ing into a stable fixed point, as rj commutes with -ffpAM 
in the limit £ — -> 00. As a consequence the differentials of 
the one-particle energies vanish in this limit, which can 
be seen on the basis of Eqs. ((TC|l - lf^ . 

The numerics generally comprises two tasks, one of 
them being the numerical integration of the flow equa- 
tions, which is done by a Runge-Kutta integration of fifth 
order with adaptive step size control. The second part of 
the numerical calculation is the self-consistent evaluation 
of the expectation values {Pf), {n^„) and (Ako-). 

Since the differential equations are coupled and possess 
a dependence on k-dependent expectation values, the nu- 
merical integration turns out to be a rather arduous task. 
Consequently, we restrict ourselves to the solution of a 
one-dimensional system. As we consider merely the para- 
magnetic phase of the PAM, we believe that one dimen- 
sion is sufficient to examine the most general aspects of 
heavy-fermion behaviour. 



V. RESULTS 

In the following section we present the results of the nu- 
merical calculations for the PAM. In a first step we anal- 
yse the flow of the one-particle energies of the Hamilto- 
nian. All calculations are performed for systems of 1000 
sites in one dimension. As we have mentioned above, we 
merely consider a dispersion-less hybridisation. Notwith- 
standing, the hybridisation Vk shows distinct flows for 
different wave vectors k due to its dependence on the 
electron energies (cf. Eq. l\W\i ). For the conduction elec- 
trons we use a constant density of states (DOS), i. e. 
a linear dispersion in one dimension, and we consider 
T = 0. 

For the upcoming discussions we use lattice constant 
of a = 1 and W ^ 2 for the bandwidth of the bare 
conduction band, so that all energies are given in units of 
the half band width, which is usually of order 1 — 10 eV ~ 
lO"' — 10^ K. As we consider one-dimensional systems, the 
wave vector is rather referred to as wave number. 



8 




FIG. 1: Flow of the quantities Ek — £k for subset wave numbers 
near the critical wave number (see text) for a system with 
4 = 0.9, V' = 0.3 and Vf = 2. 
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FIG. 2: Flow of the quantities V = 1/N J^k '^k (black solid 
line, left scale), Vk for the critical wave number (black dashed 
line, left scale) and e/ (blue solid line, right scale). The system 
parameters are as in Fig. 



A. One-particle energies 

In the following discussion we address the one-particle 
energies e^, £/ and Vk. For a particular parameter regime 
their flow diagrams are depicted in Fig. QJand Fig. |2] For 
the electron energies it is more convenient to investigate 
the difference to their initial value, namely ek(^) — 6\^, 
as a function of the flow parameter i. This quantity is 
shown in Fig. ^ for a range of wave numbers near the 
crossing point. The latter is characterised as the par- 
ticular wave number, that divides the energies between 
increasing and decreasing quantities. This behaviour can 



be understood on the basis of Eq. H16(l . as the sign of 
the right hand side depends on the wave number. For 
the system in Fig. ^ we find positive energy differences 
£^ — ejj. above a critical wave number of k = 0.4997r. As 
the graph shows, all energies start off with small modifi- 
cations up to values of £ ~ 1 before they enter a region 
of significant rise or fall. In the region of £ > 10^ they 
have already converged to their limit value. The most re- 
markable result of their behaviour is definitely the sharp 
crossover between increasing and decreasing energies. 

A much simpler picture is provided by the graphs in 
Fig. 121 As demanded, the hybridisation energy Vk van- 
ishes for all wave vectors for £ > 10^. In Fig.|3the flow of 
the quantity V = l/-/V^^Vk is shown (sohd black fine) 
beside the flow of the hybridisation Vk for the critical 
wave number k = 0.4997r. For the latter the onset of a 
significant modification is delayed to higher values of £, 
and it drops sharply to zero in the region of ^ ~ 10^. This 
suggests that the Vk with wave numbers near the upper 
and lower band edge are renormalised at lower values of 
£. Altogether we can conclude that all hybridisation el- 
ements can be considered zero above £ ~ 10^ . On the 
other hand the one-particle energy which characterises 
the local / level sees a sharp increase in the region of 
£ ^ 10° — 10^ and converges rather fast. The numeri- 
cal results show that the present choice for 77k triggers a 
smooth and continuous evolution of the one-particle en- 
ergies of the Hamiltonian 113|) . The convergence to their 
limit values occurs on a characteristic scale in the region 
oi £ > 10^. For the present results we use the bound 
V = 10^ ^'^ as convergence criteria for the numerical inte- 
gration of the flow equations. This corresponds to a value 
of the flow parameter £ ~ 3000 — 5000 in the considered 
parameter regimes. However, for a better clearness the 
integration for the system in Fig. |21 is carried out up to 
a value oi £ — 10'*. 

On the one hand we have a rather good convergence 
of the flow equations due to their smoothness and their 
evolution to a stable fixed point. On the other hand the 
self-consistent evaluation of the expectation values is re- 
stricted to certain parameter regimes. From the SBMF 
theory we know that it cannot describe the phase where 
all / sites are singly occupied, i. e. (P°) = 0. Like- 
wise, in the present approach we only obtain valid results 
for a finite number of empty / sites, and the solution 
breaks down as (Pj^) vanishes. This breakdown coincides 
with the energy e*f dropping below the chemical poten- 
tial. However, it is possible to obtain solutions near this 
critical point. The limit of singly occupied / sites can be 
approached either by decreasing the hybridisation energy 
or by increasing the difference ^ — e\. With a fixed value 
for £\ the solution breaks down below a critical value of 
V^ , and, similarly, this point is reached for a fixed V^ 
below a critical £\. Stable solutions for small expecta- 
tion values (P°) can be obtained, preferably, at rather 
moderate values of /i — £\ and small V. This restriction 
is believed to be lifted up by considering corrections of 
order i^7^ 1251. 



Within the above factorisation scheme it is not suffi- 
cient to incorporate the lyj contributions simply by con- 
sidering them in the Eqs. (|16|I - H2UI) . Indeed, it turns out 
that this procedure would lead to unphysical solutions 
escpecially in the case of small i^/. On the other hand 
the dynamics of the / electrons is even more affected 
by these corrections and cannot be properly calculated 
in the limit of small degeneracy. A consistent consider- 
ation of the Vf therefore needs an enhancement of the 
factorisation process. As we pointed out in Sec. IIII Al 
the decoupling of the higher order interactions preserves 
their character with respect to the variation of the par- 
ticle numbers for each subsystem. This assumption is 
too restrictive in the case of small degeneracy, as charge 
fluctuations which are induced by the flow of the Hamil- 
tonian become more important. Obviously, the extension 
of the factorisation scheme in this sense would lead to a 
more complex system of differential equations and is not 
considered in the present work. 



B. Electronic properties 

In this section we present the numerical results for 
the electronic structure of the PAM. The dispersion re- 
lations of the quasiparticle energies are shown in Fig. |3| 
and Fig. |3| 




FIG. 3: Dispersion of the quasiparticles for the system with 
parameters as in Fig.0 The energies e^ and £j + (PP)AJJ. are 
represented by solid black and blue lines, respectively, whereas 
the the analytical dispersions are represented by staggered 



lines. The dotted line stands for e*f. 



The general form of the one-particle energies is quite 
similar to the results of the SBMF theory. In both ap- 
proaches the spectrum splits in two bands which are di- 
vided by a finite gap. Considering the chemical potential 
lying in one of the bands, the system possesses a metal- 
lic ground state. The situation of an insulator is also 
conceivable, if fi is right in the gap. In the following con- 
siderations we only focus on the metallic case, though the 



insulating case can also be examined within the frame- 
work of flow equations [2g . 

The quasiparticle energies in Fig.|31are obtained for a 
system with e^r = 0.9, V^ = 0.3 and Vf = 2. The chemical 
potential is at /x = 1 and the bandwidth of the conduc- 
tion electron band is H^ = 2. Since we have neglected 
terms of order vj within our calculations, the presented 
results for a twofold degenerate / level are deemed an 
extrapolation down to small values of Vf. This kind of 
assumption has been frequently used in terms of various 
l/vf expansions in order to describe physical relevant 
systems. The excitation energies ej^ (solid black line) and 
£*f -\- (P°) (solid blue line) form two continuous bands, 
that are separated by an indirect gap. The latter deno- 
tation refers to the fact, that the minimal excitation en- 
ergy from the lower to the upper band occurs at different 
wave vectors. The renormalised energy e*f (dotted line) 
lies right in the centre of this gap. The self-consistently 
evaluated results for this energy and for the expectation 



value of the empty / sites are obtained as e*i 



1.032 



and {Pf) = 0.285, which states, that the system is well 
within the mixed-valent regime. Though the results from 
the SBMF theory of e} gg = 1.111 and (if )sb = 0.577 
also suggest a mixed-valent system, the latter is signifi- 
cantly further away from the integral-valent state. This 
tendency is a general result of the flow equations ap- 
proach and can be better understood in the context of 
a characteristic energy scale (cf. Sec. IVCl belowl. The 
quasiparticle energies obtained from the numerical solu- 
tion can be compared to those of the analytical treatment 
in Sec. IIV Al Their analytical expressions are given by 
Eq. (|38|l and Eq. H40|l. which have proved to be formally 
equivalent to the corresponding results of the SBMF the- 
ory. For the expectation value {Pf) and the energy £% 
we have used the numerically evaluated results in Fig. |3| 
The graphs show a qualitatively good coincidence of both 
results for the excitation energies. However, near the 
crossing point of c-like and /-like excitations the devi- 
ations become most significant. Especially, the upper 
band shows a slight depletion near the crossing point, 
which results in a Van Hove-like singularity in the quasi- 
particle DOS. This feature is a result of the consideration 
of a renormalisation of the energy Ef on a large £-scale, 
which is driven by local fluctuations. The analytical so- 
lution is based on the assumption of a fast flow of this 
energy to its limit value, which is contrasted by Fig. ^ 
and Fig. [5] Here we perceive the characteristic ^-scale 
of all energies in the same region. However, the quasi- 
particle energies of the analytical solution show a good 
correspondency to their numerical counterparts. 

Whereas the system in Fig. O was classified as mixed- 
valent system, Fig.^jis very close to the integral-valent or 
Kondo behaviour. This manifests itself in the very small 
expectation value {Pf) = 0.083 and the energy e!- — 
1.003 lying only slightly above the chemical potential. 
This system also features an evident decrease of the gap 
width. In comparison with that result the SBMF theory 
delivers {Pf)sB — 0.479 and e) sb = 1-036 indicating a 



10 




FIG. 4: Dispersion of the quasiparticles as in Fig. I^for a 
system with e^f = 0.9, V" = 0.2 and Vj = 2. 



mixed- valent behaviour. This distinctive tendency of the 
flow equation result towards integral-valence is discussed 
in connection with the Kondo temperature in Sec. IV CI 
below. 

The drastic change in the character of the quasiparti- 
cles from c to / character is a key result of the numerical 
calculations, as this contributions evolve from the bare 
conduction electron band and the dispersion-less / level, 
respectively, during the flow. This fact also ensures the 
right result in the limit of a system without hybridisa- 
tion. Notwithstanding, the different types of excitations 
form two continuous bands. In the context of the SBMF 
theory the side conditions are usually chosen in the way 
that each quasiparticle bands carries a single character. 
The switch of the excitation character in each band is 
also reflected by a change in the behaviour of the weight- 
ing factors a\^ and /^j^. The dispersion of their squares 
is shown in Fig. [S] As the evaluation of the expectation 




FIG. 5: Weighting factors (a^)^ (black sohd line) and 
{P?){l3ti)^ for the same system as in Fig. |3] 

values in Sec. IIIIBI incorporates both the quasiparticle 



energies and the weighting factors, these quantities are 
not affected by the switch of the excitation character. 



C. Kondo temperature 

Historically, the Kondo temperature Tk was intro- 
duced as a low energy scale associated with the screen- 
ing of a magnetic impurity in a metallic host. The mi- 
croscopic models for the description of such a situation 
are the single impurity Anderson model (SIAM) and the 
single impurity Kondo model (KLM). Renormalisation 
group calculations H^ as well as SBMF theory 0, H^ 
show a dependence of this energy scale on the Kondo 
couphng J = (T^')^/|eV — ^| of the form 



fceTK oc exp 



JvfPo 



(44) 



This relation states, that this energy scale vanishes expo- 
nentially as V^ goes to zero, a result that cannot be ob- 
tained by ordinary perturbation theory in V^ . In a peri- 
odic array of magnetic impurities beside the conventional 
screening of the local moments a coherence scale T* seem 
to exist, which is referred to in the literature as coher- 
ence temperature, indicating the onset of Fermi-liquid 
behaviour. Whereas the Gutzwiller calculations for the 
PAM suggest a significant enhancement of T* compared 
to the Kondo temperature of the SIAM, the SBMF the- 
ory provides a coincidence of both temperatures for small 
T* . The comparison of both energy scales has been ded- 
icated a controversial discussion. While some numerical 
works show an enhancement of T* for systems at half 
filling 11] , a more qualitative discussion by Nozieres |23| 
relates this issue to the so-called exhaustion problem. As 
Nozieres points out, only a fraction of riser = Po^k con- 
duction electrons contributes to the screening of local 
moments. As a consequence of this, he deduces an up- 
per bound of T* ~ p^T^/nf, where n/ is the number 
of magnetic impurities, for the lattice energy scale. This 
behaviour is supported by numerical calculations by Vid- 
hyadhiraja et al. |l,l |. 

In accordance to the SBMF theory [2^ we define the 
Kondo temperature by the difference of the renormalised 
/ level and the chemical potential as 



ksT* ^ e*f - p 



(45) 



As we have mentioned at the beginning of this section, 
our method provides valid results in a certain parameter 
regime. By lowering the hybridisation y or by increasing 
the difference p — e^^ the self-consistency scheme breaks 
down as the expectation value {P^) vanishes. The van- 
ishing point coincides with e^ slipping below the chemi- 
cal potential, which is consistent with the assumption of 
positive results for the lattice Kondo temperature. 

The evaluation of the Kondo temperature also allows 
us to estimate the resulting effective masses of the quasi- 
particles by comparing this energy scale to the band 
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width of the bare conduction electrons. This leaves us 
with the relation m* /m ~ W/k-^T*. If we consider the 
system in Fig. 13 we obtain an effective quasiparticle 
mass of m* /m ~ 60, indicating the system being in the 
mixed- valence regime. The enhancement of the quasi- 
particle mass is due to the flatness of the quasiparticle 
band near the Fermi surface, and it is also associated 
with a large quasiparticle DOS. For the integral-valence 
system in Fig. ^ we even deduce an effective mass of 
TO*/m-600. 

To address the issue of the functional trend of the 
Kondo temperature for the numerical solution, we evalu- 
ate this quantity by varying the coupling J. In particular, 
we examine the dependence of T* on the hybridisation 
V\ The result for the system with e\. = 0.9 is shown in 
Fig. El Here V^ is gradually decreased until the solution 
breaks down below the value V^ = 0.11. The numerical 
results are represented by circles, and the blue solid line 
characterises the fit f3exp{—S/Jvfpo). Our findings for 
the fit parameters are f3 ~ 0.20 and 6 — 1.65. The latter 



conduction electrons this yields 
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(46) 



Likewise we arrive at a representation for the /-electron 
Green's function: 
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The electronic structure can be made descriptive by eval- 
uating the DOS of the c and / electrons. These quantities 
are given by p^uj) = -1/{t:N) J2)^ lm{{c^^; cl^)) {uj + iS) 
and similarly for the / DOS. The results are shown in 
Fig. |7|for the same system as in Fig. |3| Here we have 
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FIG. 6: Kondo Temperature as a function of 1/J (see text) 
for a system with £j = 0.9 and Vf — 2. The blue solid line 
represents the fit /3 exp{~S / Juf po) with /3 = 0.20 and 5=1.65. 

identifies the value of the exponent and is significantly 
enhanced to the corresponding value for the SIAM, i. e. 
5 = 1. That means that the lattice energy scale features 
a remarkable reduction. As we have mentioned above, 
this behaviour gives credence to the exhaustion scenario. 



D. Green's functions 

In Sec. lIIBl we have explained the evaluation of correla- 
tion functions within the framework of the flow equation 
method. As such quantities are invariant with respect to 
the unitary transformation, we are in the advantageous 
position of describing thein with respect to a simple ef- 
fective Hamiltonian. In this context we use the results 
for the transformed operators from Sec. IIII Bl For the 



FIG. 7: Densities of state pc{i^) (black solid line) and Pf{uj) 
(blue solid line) for the same system as in Fig. |3 

used Lorentzian broadening with a width of s = 10""*. 
The DOS of the c electrons remains mainly flat as we 
expect from the constant DOS for the bare conduction 
electrons, but it shows a gap around the renormalised en- 
ergy e*f. On the other hand, the DOS of the / electrons 
increases sharply around the gap, meaning that there is 
a high DOS at the Fermi energy {p = 1). This picture 
qualitatively coincides with the electronic structure ob- 
tained from other approaches to the PAM, such as the 
SBMF theory, and conveys the idea of heavy quasipar- 
ticles. A remarkable difference to the SBMF result is 
an additional excitation peak in the upper band, which 
is due to the above discussed deformation of the quasi- 
particle bands (cf. Fig. 13). It possesses the character of 
a Van Hove-like singularity, as the DOS becomes very 
small in this energy range. This large DOS comes from 
the consideration of local fluctuations, which lead to a 
renormalisation of the / level. This can be contrasted to 
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the case, in which we keep this energy constant during 
the flow process. With this assumption the quasiparticle 
bands coincide with those of the SBMF theory and do 
not show this type of excitation. 



VI. CONCLUSIONS 

In this paper we have presented the description of 
heavy fermion behaviour within the framework of Weg- 
ner's flow equation method. This work provides a new 
semi-analytical approach to HFS. The basic idea of this 
scheme is the derivation of an effective Hamiltonian, that 
describes the essential physics of the system. Though this 
Hamiltonian is not diagonal in terms of usual fermion 
creation and destruction operators, it characterises two 
decoupled subsystems of electrons. Due to the Hubbard 
operators we have used for the / electrons, the according 
subsystem needs further approximations for the evalu- 
ation of physical properties, such as expectation values 
and correlation functions. In general, our findings con- 
firm the picture that is provided by previous calculations, 
such as the SBMF theory. However, the flow equation 
approach to the PAM yields an additional excitation in 
the upper quasiparticle band and a significant decrease 
of the Kondo temperature. 

In order to obtain a closed set of differential equations, 
we had to resort to a decoupling approximation, and we 
have only considered contributions in the leading order 
of Vf . The idea of the second step is similar to the proce- 
dure used in the SBMF theory. However, in this context 
the two approaches cannot be compared in a one-to-one 
fashion. As we have seen in the discussion of Sec. IV Al 
an improvement of the present approximations can be 
obtained by an enhancement of the factorisation scheme 
due to the importance of charge fluctuations in the case 
of small degeneracy. An important difference of the flow 
equation approach to the SBMF theory is the determina- 
tion of the energy £*s. Whereas in the SBMF theory this 
energy is obtained by a mean field procedure and has 
to be determined self-consistently, the present method 
uses the integration of a differential equation. The proce- 
dure of self-consistency within the flow equation method 
only refers to expectation values. As we have shown in 
this paper, the differential equations can be solved both 
analytically and numerically. Though the flrst was per- 
formed on the basis of further simplifications, and it was 
intended to give a rough estimate, it turned out to coin- 
cide qualitatively good with the numerical solution. In 
accordance to the SBMF theory the present calculations 
show a restriction to parameter regimes where the oc- 
cupation number of the /-sites is less than unity, as the 
solution breaks down in the limit {Pf) — 0. However, we 
obtain good results for quite small deviations from this 
point. As we have stressed in Sec. lVAl this breakdown is 
a result of the scheme of selfconsistency, while the inte- 
gration of the flow equations is rather robust and shows 
a good convergence even in this limit. 



In the present paper we studied the paramagnetic 
phase of the one-dimensional system. From recent works 
we know that the ground state of the KLM in one di- 
mension is either ferromagnetic or non-ordered. On the 
other hand DMRG studies for the C/ = oo PAM also 
showed regions of short-range antiferromagnetism in the 
Kondo regime [I3, whereas the ground state is found to 
be paramagnetic for the mixed-valent region. For the lat- 
ter our results are consistent to those obtained in Ref.llTl 
Near the Kondo regime magnetic interactions, such as the 
RKKY interaction, become more important. For a bet- 
ter understanding these contributions must be incorpo- 
rated in the present scheme. Due to the semi-analytical 
character of the flow equation method we believe that 
this task could be accomplished for larger system sizes 
as usually considered in the DMRG. Another intriguing 
issue which was addressed in recent works about one- 
dimesional heavy fermion systems is the study of Lut- 
tinger liquid behavior 29]. However, this subject was 
merely discussed for the KLM and an occurence of a Lut- 
tinger liquid in the PAM is still an outstanding problem. 
The examination of this issue though needs a further en- 
hancement of the calculations and cannot be obtained on 
a simple mean-field level. 

Our findings for the excitation of the quasiparticles 
allow a clear distinction of c and / contributions. Gener- 
ally, the electronic structure of the PAM comprises two 
quasiparticle bands, which are divided by a gap. The 
DOS is strongly enhanced at the Fermi level, resulting in 
quasiparticle masses up to several hundreds of bar elec- 
tron masses. Such a behaviour characterises the exper- 
iments for heavy fermion compounds. In comparison to 
previous theoretical works on the PAM, the present pa- 
per shows an additional Van Hove-like excitation peak in 
the upper quasiparticle band, which is due to a significant 
fiattening of the dispersion. 

We also have investigated the dependence of the Kondo 
temperature on the hybridisation V^ . This quantity van- 
ishes exponentially for small values of V^. Our findings 
for T* are consistent with the qualitative arguments pro- 
posed by Nozieres, suggesting a significant decrease of T* 
compared to the Kondo temperature of the SI AM. 
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APPENDIX A: DECOUPLING SCHEME 

Here we present the factorisation procedure discussed 
in section IlII Al for the contributions in Eq. (|14|l . With 



the denotations {c\^„c-^^„) 



^L) and {fl,U 



we obtain 



^EE^'^'"'^"'^>^^-4.^.'.^°+h.c. 



kk' i(T 



kcr 



kcr^kcr 



A^ 



E^k^kKJi^° 



n ("L> 



(Al) 



kCT 



^EE^^^''"'^"'^'^^-4.'^k'.'/l'/..+h. 



kk' zctct' 



2E^kl4(nL 



/\.t 



cl„c. + 



■E^'kVkKj/lA 



kfj 



ik(T 



2E'7kli(nt)K,) 



(A2) 



kCT 



In order to obtain the appropriate contributions to the 
one-particle terms of the Hamiltonian from Eq. lAll it is 
convenient to make use of the relation P° + ^^ f tafia = 
1. The remaining interactions are decoupled as follows 



kcr i^j 



J2v^.^^{P^){cUua+fl 
kcr 

1^77kAk(Ak.)/'° 



ka'-kCT; 



N ■ 



ikcr 



zkcr 

-^ E E^kt„e-''^-4,44J,,, +h. c. 

If iT iT ' 7 =;^ T 

kcr 

+ ^E^kAk(Aka)/L/^. 



(A3) 



kcrcr' i^^j 



ikcr 



E'?kAk(ni>(^ka) 



(A4) 



ka 
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where (Ak^) = (cL/k^ + /Lck^) ^as used. 



APPENDIX B: LOCAL EXPECTATION VALUES 

The expectation value (P°) is best evaluated by mak- 
ing use of the the relation 



P° = 1 - V f ^ f . 



(Bl) 



Further we can take advantage of the translations invari- 
ance of the sjstem, which allows us to write (fJafia) ~ 
1/-^ Si ifjafja)- Consequently, it is convenient to regard 
the unitary transformation of the operator 



F=l/Nj2flU 



(B2) 



Its ^-dependence can be described in the usual way by 
dF{£)/d£— [i]{£),F{£)], where the generator is given by 
((Tn|l . Thus the following ansatz for the operator is expe- 
dient 



kcr 



-kCT^kCT 



Wa 



kCT 



a'^kaj 



m 

N 



(B3) 



In this connection the boundary condition A^^'' (^ = 0) = 1 
hold, whereas all the other coefficients are equal to zero 
for £ = 0. Evaluating the commutator we arrive at 



['?' ^] = - ]^ E ^^i-B^ + ^*'^ + <^")4'^) (4cr/k. + flaC^a) + ^ E E 2^kGk(nL) 4 /.. 
ko" ik era' 

+ ^ E 2^kGk(P°)cLck. - ^ E 2^kGk/L/k. + ]^ E 2^kGk(nL> (l - (^°» 

kcr kcr kicr 



Here the decoupling approximation, which was discussed 
in Sec. IIII Al is already done. As a result we obtain a 



^E2Kc.)^/kGk 



(B4) 



I 

closed set of flow equations for the coefficients 



(B5) 
(B6) 
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The following exact relations can be derived 



AW=l + lV(nL)A, 



(2) 



A 



(2) 



ko- 



{Pf) 



Furthermore the differential equation 

-dBk 



/pO\ dGk 

^ '^ d£ 



d£ 



Ad): 



d£ 







(BIO) 

(Bll) 
(B12) 

(B13) 



holds. As a simplification we regard the summation on 
the right hand side of (|B10|1 a small correction, and 
we neglect it. Such an approximation is well confirmed 
by numerical simulations. As a consequence we obtain 
^(1) — \ for the entire flow, and the differential equation 
(jB13p can be integrated straightforward 



d£ 



arcsin[l - 2Bk] - -2J{P^)rjk 



(B14) 



With the help of the unitarity relation H25() and l|26|) we 
can identify _Bk = (P°)/3^, and the expectation value 
thus reads 



(P' 






(i-Kf)*(cLO 



+ i^-m')*{fLLr 



(B15) 



An equivalent procedure can be used in order to derive 
the expectation values (cl^cr^ka) ^^'^ i^Wa)- Due to their 
non-local character the results coincide with those ob- 
tained in Sec. IIII Al Therefore it is more reasonable to 
use the single operator transformation, because this ap- 
proach can be used to describe correlation functions, as 
pointed out in Sec. IV Dl 
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